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The very small size of optical nonlinearities places strict restrictions on the types of novel physics one can 
explore. For an ensemble of multilevel systems one can synthesize a large effective optical nonlinearity using 
quantum coherence effects but such non-linearities are technically extremely challenging to demonstrate at the 
single atom level. This work describes how a single artificial multi-level Cooper Pair Box molecule, interacting 
with a superconducting microwave coplanar resonator, when suitably driven, can generate extremely large opti- 
cal nonlinearities at microwave frequencies, with no associated absorption. We describe how the giant self-Kerr 
effect can be detected by measuring the second-order correlation function and quadrature squeezing spectrum. 

PACS numbers: 03.67Hk, 03.67Lx, ()5.50.+q 



Introduction:- Given a sufficiently large optical nonlinear- 
ity with low quantum noise, it should be possible to gener- 
ate and observe strictly quantum effects in electromagnetic 
fields. Examples of such effects include quadrature squeez- 
ing H], generation of a superposition of macroscopically dis- 
tinct quantum states [2|, optical switching with single pho- 
tons |3| and measurements of nonlocal correlations of entan- 
gled photon states lH. So far, the successful demonstration of 
these effects has been limited to implementations with pho- 
tons and atoms. The main obstacle for such an implemen- 
tation -spontaneous emission- can be bypassed by exploiting 
quantum coherence effects in multilevel atoms. Such effects 
include coherent population trapping (CPT) I5j, electromag- 
netically induced transparency (EIT) |6l and others. 

Recently, a novel system was shown to be capable of im- 
plementing basic quantum optical systems. Circuit quantum 
electrodynamics (cQED) is an on-the-chip counterpart of cav- 
ity QED systems Q, that employs a quantised microwave 
mode held in a Co-Planar Waveguide resonator (CPW) (sub- 
stituting the standing-wave optical cavity) and a Cooper Pair 
Box (CPB) (instead of two-level atom trapped in the cavity). 
This system offers an unprecedented level of tunability and 
flexibility in the implementation of strong-coupling interac- 
tion limit. The aim of this paper is to show that, by designing 
and utilizing multi-level coherent processes in "artificial su- 
perconducting atoms", Circuit-QED systems can display ef- 
fects completely analogous to CPT and EIT. In this work we 
show that by designing a specific type of four level "arti- 
ficial superconducting atom", and arranging for EIT in this 
system when coupled to the CPW, one obtains a giant Kerr 
nonlinearity (Fig. [TJ. EIT is based on the use of dark reso- 
nance where quantum interference cancels the absorption of 
the probe signal. A strongly detuned fourth level provides an 
ac-Stark shift to the ground state |3), resulting in a self-Kerr 
nonlinearity free of spontaneous emission noise. Schmidt and 
Imamoglu |8| predicted that this iV-scheme can give rise to 
several orders of magnitude enhancement in Kerr nonlinear- 
ity as compared to conventional schemes. Their prediction 
has been verified in a recent experiment by Kang and Zhu 
lIPl, although in a semiclassical regime. In our work we show 



|4)~|TT> 




Figure 1: N-system fSl, constructed from two coupled two-level CPB 
systems where transitions |1) ^ |2) and |3) <-> |4) are coupled with 
strengths g\, g2, to the microwave photonic mode o of frequency 
LJa, held in the CPW resonator. The transition |2) <-> |3) is driven by 
the semiclassical control field flc, while 7ij, denotes the decay rate 
from |i) \j). We define A = los4 — uja and 5 — ul>\2 ~ ^a- 

that, due to the very strong coupling between the "artificial 
superconducting atom" and the CPW, the resulting giant Kerr 
nonlinearity is predicted to be three-to-four orders of magni- 
tude larger than what has been so far experimentally demon- 
strated. Moreover we predict that the effect could be quite 
robust against dephasing. The presence of such large Kerr 
nonlinearities in a high-finesse cavity could result in complete 
photon blockade giving an effective two-level behavior for the 
cavity mode fTOl [TT] [El [El Pi- A large effective self-Kerr 
nonlinearity allows for a photon blockade effect, i.e. a source 
of well resolved single photons. It could also be potentially 
adopted for use as a microwave single photon detector and to 
implement conditional quantum logic. 

Below we show that one way to tailor the required multi- 
level "artificial superconducting atom", giving a four-level 
A^-scheme (4 atomic levels with transitions in the shape of 
the letter A^), is via the straightforward capacitive coupling of 
two CPBs. When this system is coupled to a quantized CPW 
field, severe photon blockade can be observed, corresponding 
to a huge nonlinear photon-photon interaction. The effect of 
quadrature squeezing is also predicted. Together, the Giant- 
Kerr effect of the N-scheme, with the very strong coupling in 
cQED (due to the large dipole moment of the CPB and the 
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small mode volume of the CPW), gives the posibility of im- 
mensely large optical nonlinearities at the single photon level. 

CPB Molecule:- Normally CPBs are operated at the charge 
degeneracy point and act as an effective two level system 
when the charging energy greatly exceeds the Josephson en- 
ergy, i.e. Ec ^ Ej. To model a multilevel atomic sys- 
tem, one can capacitively couple two CPBs together to form 
a CPB molecule. For weak coupling, the CPB-molecule's 
states III), are nearly degenerate, while for large coupling 
the corresponding eigenstates are non-perturbative superposi- 
tions of these bare states and are strongly split. This results 
in a formation of multi-level system ifTSl . For zero detunings 
A = (5 = 0, we arrange for resonance between the 34, 12 tran- 
sitions and the cavity, ^^34 = UO12 — ^a- We now show that by 
operating the capacitively coupled CPBs at the co-resonance 
point, a symmetric A^— system is realised (see Fig. [TJ. With 
this design we also have the important flexibility to adjust the 
level structure: by tuning the flux threading the CPBs equally 
within the CPB-molecule, an asymmetry is introduced leading 
to nonzero detunings A. S. 

Considering two capacitively coupled CPBs (Fig.|2|A)), the 
associated Hamiltonian is a combination of "Kinetic" (KB) 
and "Potential" (PE) terms associated with the phase dif- 
ference of the wavefunctions across each Josephson junc- 
tion 01(2). From the properties of Josephson junctions 
one has (pj — 2eVj/h, where Vj is the voltage drop 
across the Josephson junction, and the current through 
the junction is Ij = IcSiiKpj. The total energy required 
to charge all the capacitors in the circuit gives KE = 

1 /2 Ej=i (C^^'^ V''^^ ^ + Cg ^ Vg'^ ^) + 1 /2C('") ^, while 
flie total PE = Y,l=i - cos (/.C^)). sf'^ is the Joseph- 

son energy of the fc* junction. Starting from Kirchoff 's laws, 
we obtain the classical Hamiltonian 
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H = Y, [Ei^\n('^ -N^gJ^f +E':p{\- cos <j)^J^)^^ 

+i?„(nW-7V«)(n(2)-7Vf ) , (1) 

where the following quantities have been defined 

{E^\E,r.} = 0.5(2e)2{(c(?,)-\(C7j;)-n, 

and where the number of excess Cooper Pairs on the gates 
is nI^^ = -C^^V(jV(2e). We consider the quantised 
version of (fill and focus on the low energy dynamics by only 
including |riJ^\ n^^^), rS^'> = 1,2 states in an expansion of 
the Hamiltonian operator Setting (5'^-* = N^f '' — 1/2, and 
denoting Z = |0)(0| - |1)(1|, X = |0)(1| + |1)(0| gives 

H/h - - + JZ(1)Z(2) , (2) 

where hwi^^ = E^^^ + £;„/2, hw'^P = s'^p /2, and hJ = 
Em/4:. If the individual CPBs are identical then the individual 
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Figure 2: (A) Circuit of two capacitively coupled CPBs with individ- 
ual gate bias; (B) Schematic of possible physical arrangement of an 
identical = ^r), pair of CPBs to yield a CPB molecule. 

Josephson energies in Q can be modulated by local flux tun- 
ing via hwi^^ ^ E^p = 2£jcos(7r$(jV*o)- Working at the 
CO- resonance point 6^^^ = 0, and setting ujP' = lo'P^ = w^; 
yields 



(3) 



The eigenenergies of ([3]l, labeled as in Fig. [T] are 
{Ei,E2,E3,Ei) = J X (e,l,-l,-e), where = 1 + 
Auj^/,P. The transition frequencies are then LO21 = LO42 = 
J(e + 1) and W23 = 2J. To arrange for nonzero detunings 
A 7^ (5 (i.e., R = 00^2/^21 7^ 1), one must move slightly off 
the co-resonance point by introducing equal strength Zeeman 
terms in (|2|, 
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with uj'i^^ — Lo'^P 5^^\ R can be set to any value i? > 1 so 
adjusting and lo^ fixes the detunings A, 5 to any desired 
value without individually addressing each CPB. 

To estimate the size of the coupling J, consider a closely 
spaced pair of CPBs as shown in Fig.|2];B) il6i . The capac- 
itance between the two CPBs can be expressed as C'™) = 
27rreeow/e where ^^K{^) + K{^)- K{^^^), with 
K{x) — a; sinh~^(l/a;) + sinh^^(x) and e = 9 the relative 
permittivity of the substrate material fVl]. We can estimate 
that if {I, w, r) = (50, 10, 0.5)/im we obtain a large splitting: 
{ujy\ J) = (14.7, 2.3)GHz. Due to the nonlinear relation- 
ship between C'™) and hJ = Em/ 4, the frequency W32 ^ 2 J 

for the co-resonance case can be set to be MHz-GHz depend- 
ed ^\ 2") 

ing on the CPB geometry. The value Em'wx — 2ujx ^ ^ 8 
GHz |7 1 can be reduced via adjusting the flux threading both 
CPB (assuming identical CPBs = ^r)- As a prototype. 
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Table I: Table of ri from Eqn |5j, measuring the size of the effective self-Kerr nonlinearity for various quantum systems coupled to cavities. We 
take the quoted values for [g, k, 7), and consider additional driving (see Fig. [T|( such that {g/Q,c)'^ = 0.1, [condition of validity for Eqn (|5|], 
take (A, 5) — (7, 0) and 7; = 7. The effective self-Kerr nonlinearity for the CPB-molecule/CPW system has the potential to be extremely 
large. The actual model studied in this Letter contains extra decays in addition to those traditionally studied in A'^— systems. 



we choose UJ43 = UJ21 = o^a = 5 GHz, giving vanishing 
detunings A = (5 = 0, if one arranges to work at the co- 
resonance point = 0. Then by driving the 2 <-> 3 tran- 
sition with RF at frequency 0-123, one gets J = 0.2GHz, and 

(1 2) (1 2) 

ujx — Ej =2.6 GHz. To arrange for non- vanishing 
detunings for example, we can choose to set Wa;/27r — 4GHz, 
J/27r = 0.2GHz, = Jf^ = 27r x 16GHz, and with 
«;/27r = lOOkHz, one can obtain significant differences be- 
tween the transition frequencies R = (JJ34 — iO\i, to yield 
Rj K ^ 1. One can achieve these moderate detunings even 
when the gate bias is offset from the precise co-resonance 
point by a tiny amount (5^^^ = (5^^^ = 2.8 x 10^^. In the 
special case of identical CPBs operating exactly at the co- 
resonance point one has selection rules which could alter the 
transition strengths within our system but as we work off co- 
resonance we expect all transitions and decay paths to be al- 
lowed. This is a significant departure from the atomic A^- 
system studied in |8|, and we will include all these decay 
routes in our full model below. 

Giant Kerr Nonlinearity:- It was shown ifTOl fTTl [121 [T3l 
that an A^-system similar to one described in Fig. [T| yields 
an effective Kerr nonlinearity. Due to the absence of selec- 
tion rules, we consider all possible decay paths with rates 7^ . 
Drummond and Walls 1 18 1 analysed the pure x^'^-' (third-order 
optical nonlinear susceptibility [IJ, Sect: 5.4), Hamiltonian 
iieff = hrjo^^a? in the presence of a dissipative cavity and 
using this we will estimate the size of the effective self-Kerr 
nonlinearity r] we observe in our effective A^-system. For 
large 77, the system emits photons in an antibunched manner, 
with large waiting times between single photon emissions. 
To probe these effects, the second order-correlation function 



a) {t)a} {t + T)a{t + T)a{t)) / {a) {t)a{t))\\ will be 



calculated, in particular (/'^^^(O), for the weakly pumped cav- 
ity. This quantity is accepted as a good measure for photon 
blockade lfT0l[T2l[T4l[T9l . and allows for a direct comparison 
with the analytical expression obtained in 1 18 1. 

Before examining the combined resonator-CPB molecule 
in detail we can find a rough estimate for the achievable non- 
linearity rj using parameters from a recent optical turnstile ex- 
periments (see Table ), taking IJ2j 
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V7l3 + A^ (721 + 723)2 + (52 
which holds in the limit of {gi/VtcY < 1 ifTTl . 



(5) 



For the resonator-CPB molecule system we set 721 — 
723 = 743 = 7 = 1/^1' where Ti is the lifetime of the 
single CPB excited state, and we take Ti — 10/xs, n — IMHz, 
and choose (A, (5, fie) = (10,0, 1)7. We choose the cross- 
relaxation rates 742;3i;4i to vanish (only for the purpose of 
this estimate) and gi — §2 = g = 300k, to obtain the enor- 
mously large Giant self-Kerr strength 77 w 10 x IO^k. Al- 
though (jsjl only holds when gj/^c ^ 1. this computation 
gives one the hint that the system (with all decay channels op- 
erational) may yield very large nonlinearities and in the next 
section we model the full system to numerically verify this. 

Numerical Investigations:- Estimating the size of the effec- 
tive nonlinearity r/ when gj > is not straightforward as the 
adiabatic approximation is not permitted ifTOl . For that reason, 
we solve the master equation for the density operator p, given 
hy p— Psys + Cp, where 

Psys = -iA21 [ct22,P]-«A3i [0-33 , p] - « A41 [(J44, p] 

-igi [aVi2 + cr2ia, p] ~ ig2 [0^34 + 0-430, p] 
-i [f2*cr32 + cr23fic, p] - iEp [a + a'^ , p] , (6a) 



and Cp — tiD[a\p + 
2ABAt - {A^A,B). 
operators X'^^i) = a^j 



l^j'D[crj^]p, where V[A]B = 
Here, aij describes atomic transition 
(7ji for i ^ j, while a-jj models de- 



phasing. Ep ^ ^yPn/iUa is the amplitude of an electric field 
driving the resonator mode, where P is the power of the field 
incident on the resonator. Summation over (ij) includes all 
decay channels shown in Fig. [T] We solve p — numeri- 
cally to obtain g'-^^O). Fig.blshows the dependence on pump 
strength and (classical) coupling field strength. The coiTela- 
tion function increases with pump strength |14|. The depen- 
dence on the coupling field shows a decrease to a minimum 
(solid black line) followed by an increase. This is a conse- 
quence of a nonvanishing decay 731. Increasing fie pumps 
the population in |1) and diminishes the effect of decay 731, 
hence the decrease in g'^^^O). Subsequent increase with ilc is 
due to a reduction in the nonlinearity (j5]l. Using the analytical 
result for g^^^{0) as a function of Ep/n in [18,1, the effective 
self-Kerr nonlinearity rj/n can be deduced (shown in the in- 
set of Fig. 3). For weak driving, effective rj/K ^ — 10^ 
can be obtained. As highlighted in the introduction, using this 
"artificial A^-system", optical nonlinearities of unprecedented 
strengths are achievable using circuit QED. This result is ro- 
bust against dephasing. When we include identical rates of 
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Figure 3: Graph of log^pfg*^' (0)], for g = 300, 723 = 741 = -01, 
and 721 = 743 = 731 = 743 = .1, in units of and (A, 5) = 
(.5, .5). For low Ep and moderately large pumping Q.c, the auto- 
correlation drops to extremely low values ~ 10~®. The solid line 
indicates minimum (7*^'(0). Inset shows the effective nonlinear co- 
efficient rj vs. Ep deduced from 1181 . for the minimum 3^^' (0) (see 
text for details). The average photon number is n ~ 10"^ for maxi- 
mal value of the nonlinearity. 

dephasing on all levels (taking ^^k = Iph ^ 2.57^, as seen 
in experiments when working off the critical point), we find 
that rj is an order of magnitude smaller than what is shown 
in Fig. [3] Decreasing Ep/n further allows one to reach val- 
ues shown in the Figure (and higher), and the use of trans- 
mons 121] could reduce the dephasing rate when working off 
the critical point. 

Self-Kerr systems can produce some degree quadrature 
squeezing [22]. In Fig. [4] the spectrum of maximum squeez- 
ing (belonging to the amplitude quadrature) is shown. The 
regime where maximum squeezing is obtained [Fig. 4(B)] 
is closer to the four-wave mixing regime. Given the pres- 
ence of six coherence-destroying spontaneous decay channels, 
the result of Fig. 4 is remarkable, since pure dispersive Kerr 
medium can only achieve a squeezing of | II22I . 

While the measurement of squeezing is easy to implement 
using homodyne techniques, the measurement of second- 
order correlations involves efficient photon counters 1191 . 
which are not currently available in the MW. We note, how- 
ever, a recently technique ll23l to obtainin g'-^-'(O) by homo- 
dyne methods. 

In conclusion, we have shown that an "artificial" multilevel 
system in circuit-QED, of two capacitively coupled CPBs 
coupled to the quantized field in a CPW, leads to the gen- 
eration of giant self-Kerr nonlinearities as demonstrated by 
strong antibunching and quadrature squeezing. The size of 
the effective nonlinearity is orders of magnitude larger than 
anything previously known, owing to the combined effects of 
quantum coherence and strong coupling. 
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Figure 4: Spectrum of squeezing S{ui/k), for {A, S, g, Ep) = 
(5.13,-4.89,300,0.14) and (A): Qc = 50, or (B): ftc = 1200, 
and where 723 — 741 = -01 and 721 = 743 ~ 731 = 743 ~ -l, in 
units of K. 
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